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ABSTRACT 

Tha distinction batwaan tha usas of ralativa 
and absoiuta spactroradioaiatric calibration of ra- 
mota sansing systaBs is discussad, Tha advantagas 
of datactor*basad absoiuta calibration ara dascribad. 
and tha catagorias of ralativa and absoiuta systM 
calibrations ara listad. Tha linitationa and 
problaoa associatad with thraa coMon aatbods uaad 
for tha absoiuta calibration of rasnta sansing 
systams ara discussad. 

Tiro matbods ara proposad for tha in-flight 
absoiuta calibration of advancad sniltispactral 
linaar array (kla) systau. Ona oaJcas usa of a sun- 
illuminatad panal in front of tha sansor* tha 
radianca of which is nonitorad by a spactrally 
flat pyroalactric radicaatar« Tha othar usas a 
^*^^9** uniform, high-radianea rafaranea grotind sur- 
faca. Tha ground and aosospharic aaasuramants 
raquirad as input to a radiativa transfar pzogroi 
to pradict tha radianca laval at tha antranca pupil 
of tha ozbital sansor ara discussad, and tha 
ground instruaantation is dascribad. 

Kay vordsi Radicoatzyr Calibration! RMota Sansing. 

1. SPBCTRORADIOHBTRIC CALIBRATION— WBO w v - i™ IT? 

For purposes of this discussion, tha cooBRinity 
of reiRote sensing data users can ba divided into 
two croups I 

1. Users reouiring ralativa, but not 
necessarily absolute, spactroradicBetric sensor 
calibration. These include workers in computer- 
aided scene classification, cartographers, isiage 
processors, photointerpreters, and people concerned 
with composing large mosaics. 

2. Users reouiring absolute spectra 
radiometric calibration. These include physical 
scientists concerned with relating ground-neasured 
parameters and/or atmospheric characteristics to 
the spectral radiance at the entrance pupil of the 
space sensor. 

The distinction can be drawn that the former 
are concerned prinarilv with taking an inventory or 
comoiling a map of ground features. The latter are 
concerned primarily with understanding and 
characterizing the physical interactions taking 
place, usually through models and often with a view 
to optimizing (for example an irrigation schedule) 
or to oredicting a change in a oarticular process 


(for example an agricultural- yield) . Expressed 
differently, those concerned with relative 
calibration use data described in digital counts, 
while those needing absolute calibration use data 
referenced in terms of radiance units at the 
entrance pupil of the sensor (in this case, the 
radiometric calibration of the sensor is invoked to 
convert the digital value to spectral radiance) « 

We must emphasize that the radiance value is of no 
more worth to the average user than the digital 
value, as both refer to an incident radiance level 
at the sensor over broad, unequal spectral 
passbands, over which the spectral reflectance of 
the feature being observed can change by a 
significant but unknown amount. However, there are 
two reasons for converting the digital value to 
radiancei first, in multitemporal sensing, to 
account for anv docum en ted changes of radiometric 
calibration with timet second, to teat or utilise 
physical models in which the ground reflectance and 
atmospheric effects are measured and/or calculated 
over Identical spectral passbands as employed by 
the space sensor. 

Perfect relative radiometric response occurs 
when the outputs from all detectors in a hand are 
equal or can be adjusted during pre^ocesalng to be 
eoual, when the Incident spectral radiance is 
constant across the sensor's field of view. (Note 
that the number of detectors in a band can be as 
few as six for the Multispectral Scanner System 
(MSS) on Landsat and as many as 18,500 on future 
MLA systems.) This condition must be met 
independenUy of the spectral content of the scene. 
When this condition is not net, the image appears 
striped. If it is not scene-dependent, striping 
often can be completely removed by the histogram 
equalization methodr thus, relative radiometric 
precision can be high even thoooh the accuracy 
involved may be low. In this equalisation 
procedure the histogram of each detector output is 
compared with that of every other detector, after a 
large number of data samples (^2 X 10$) have been 
recorded. It is assumed that, if the scene is 
spatially and spectrally random, the histograms for 
a large number of samples will be identical. If 
the histograms are not identical, adjustsients are 
made during the preprocessing step to make them so. 
This procedure can be repeated for scenes of. 
different average radiance, and the relative 
responses can then be equalized over the dynMle 
range of the detectors. This procedure does not 
work if the striping is scene d e pe ndent (Ref. 1). 
Fortunately the design and the spectral band 
location for the Thematic Mapper and future MLA 
sensors will not give rise to a scene- d ependent 
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strioina nroblera of the ffaonitude of that in the 
MSS. 

Inac’eouatelv corrected relative detector-to- 
detector response causes unsightly stripino in the 
imagery. Striping can cause inaccuracies in * 
automated scene classification. In addition, even 
though the human ohotointerpreter is usually more 
foroivino than a computer, if the striping is 
severe it can cause errors in photo interpretation, 
ancorrected strioing can be aagravated by the 
apolication of some imaoe orocessing algorithms. 

Thus, it is imoortant to reduce striping as much as 
possible for image processing and classification 
putooses. The relative response of the individual 
detectors is also important for ratloing purposes. 
However, once the in-band tolerance has been met, 
the concern is with the relative stability of 
response because many applications utilize the 
comoarison of multi temporal band ratios. 

The utilization or verification of physical 
models usually reouires the use of data calibrated 
in an absolute sense. Until recently the highest 
in-orbit absolute radiometric accuracy has been 
little better than 10%. This low accuracy has been 
due to: (a) the fact that the calibration in orbit 

has often been for the focal plane onlv, not for 
the comolete system: (b) the loss In accuracy 
accompanying the transfer of calibration from the 
standard source at the national laboratory to the 
factory or . laboratory calibration site: (c) the use 
of source-based calibration procedures. 

As described later# the use of detector-based 
calibration pronises to reduce the approximately 
10% error to 1%. To what extent we can afford to 
relax or strive to exceed the 1% potential can best 
he determined from a set of well-coordinated 
measurement and modeling exercises conducted as 
part of an experimental MIA mission. The results 
of work to date on this subject are in 
disaoreement. On one hand we know that natural or 
scene variability is seldcn less than 1% and we 
also know from atmospheric modeling work (not yet 
experimentally verified) that we can expect 
considerable a ttBospheric-scatter- induced spectral 
radiance crosstalk between neighboring ground 
samples owino to the so-called adjacency or 
boundary effect (Refs. 2-5) . On the other hand, 
the sensor simulations made to date (Ref. 6) 
indicate there is an improvement in classification 
accuracy in going from 6-bit to 8-bit quantization. 
Bight-bit quantization is available on most future 
sensors. A purist miaht contend that, If 8-btt 
radiometric resolution is really justified, it 
should, for maximum utility, be associated with 
comme nsurate (<0.4%) error in absolute radiometric 
accuracy I 

In summary, we know that soectroradiometrlc 
calibration is important and that the need for 
relative or absolute calibration depends on the 
application for the data. However, we do not know 
how accurate the calibrations should be, 
particularly the absolute calibration. The moat 
encouracing aspect of the situation is that we now 
have the potential of utilizing 8-bit quantization 
resolution with commensurate absolute radiometric 
accuracy in an orbital imager at least over the 400 
to 900 nm range. This potential should be 
exploited as fullv as possible as an important 
component of any system and applications research 
aaaociated with an experimental land-observing 
system having an 8- to 10-bit quantization 
capability. 


2. OETECTOR-HAHED RADIOMETRIC CALIBRATION 

The recent work at the US National Bureau of 
Standards on self-calibrated photodiodes is 
described only briefly here. For more details the 
interested reader is referred to references 7-12. 

The calibration of the photodiode is 
accomplished by either of two differently applied 
biasinq procedures, depending on the wavelength 
region of interest. At short wavelengths, a 
negative bias is applied to remove the 
recombination centers at the Sl-SlO interface at 
the front of the detector. To do this, a contact 
Is made with the front surface using an electrode 
immersed in a conducting liquid, or the surface can 
he exposed to a corona discharge. For long 
wavelengths, a back-bias is applied to extend the 
depletion region to a depth beyond which incident 
flux penetrates. The experimental procedure is to 
irradiate the detector with a constant 
monochromatic flux level and to increase the bias 
voltage until further increase no longer gives rise 
to an increase in output signal. For both the 
short and the long wavelength ranges, the internal 
quantum efficiency saturates at a value extremely 
cloae to unity, as shown in Fig. 1. Thus the 
maximum increase in signal output obtained as a 
result of biasing can be used to determine the 
internal quantum efficiency of the detector without 
biasing, as it will be used in practice. 



Fig. 1. Typical photodiode Internal ouantum 
efficiency without biasing (lower curve) and with 
biasing (upper curve), reference 9. 

The only significant loss in the external 
quantum efficiency of the photodiode is caused by 
reflection. This can be reduced to insignificance 
by making use of three photodiodes according to the 
geoBietry sketched in Fig. 2. 



Fig. 2. A three-diode arrangement to 
minimize specular reflection losses. 


The output signals from the three diodes are 
summed to provide the signal corresponding to a 
black detector of overall quantum efficiency that 
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can be assumed to be uni tv. The second diode 
reflects the specular reflection from diode 1 to 
diode 3, which acts as a retroreflector. The 
incident flux thereby underooes five reflections, 
and at a 10% value the final specular reflectance 
is dovm to 10*^ of the Initial Incident radiance. It 
is claimed (Ref. 12) that the diffuse reflectance 
losses for clean detectors are tyoically much less 
than 1%. The photodiodes used In similar 
calibration facilities at MBS and at the University 
of Arizona, Ootical Sciences Center, are EGG UV 
444B. 

The discussion at the end of this paper refers 
to the use of self-calibrating NBS detectors for 
the spectroradlometric calibration of an MLA 
system. However, as shown in Rig. 1, the unbiased 
cuantura efficiency is wavelength dependent* 

Because of possible changes in the passband 
position of soectral filters during long duration 
space flights, it may be advisable to use the NBS 
detectors to calibrate a spectrally flat 
cyroelectric detector, at the 0. 1-0.2% level, and 
for that to be used in the in-orbit calibration. 


3. BROAD CIASSIPICATION OR CALIBRATION PROCEDURES 

Some of the moat cammonly used procedures for 
the calibration of remote sensing systosa are 
referred to In Pig. 3. The two major divisions in 
the fioure are between relative and absolute 
calibration and between the static macro-image 
response and the dynamic micro- imaae response of 
the system. Sene aspects of relative calibration 
were discussed in the first part of this paper and 
tocether with dynamic micro-image response will not 
be discussed further, beyond remarking that the 
dynamic micro- image responae is of vital interest 
in anv pixel-by-pixel analysis of remotely 
collected Imagery. The rest of this paper deals 
with the absolute calibration of remote sensing 
systems, and only the static macro- Image response 
will be considered in this context, as is usually 
the case. 

The procedures for the absolute calibration of 
a remote sensing system fall into the three 
categories shown in the bottom right of Pig. 3* 

(1) The absolute calibration of the system is 
made onlv before launch. In flight the calibration 
is checked by Irradiating the focal plane with a 


radiometricallv calibrated source and ootical 
system. The drawbacks to this orocedure are that 
any chance in the transmission of the imace-forming 
optics of the sensor system, due to the 
condensation of outoassed contaminants, will be 
undetected and the on-board calibration system Is 
also assumed to be stable through launch and 
unaffected by the vacuum, high enerev particle 
irradiation, and zero-g environment at orbital 
altitudes. The Thematic Mapper and the 
Multlspectral Scanner Svstem on Landsat-D are 
examoles of remote sensing systems calibrated in 
this manner. 

(2) The sun or an on-board calibrated source 
can be used to irradiate the focal plane through 
the image-forming optics. The drawbacks to this 
approach are the uncertainty in the knowledge of 
(a) the irradiance of the sun above the atmosphere 
and (b) the output of the calibrated source system, 
for the reasons mentioned earlier. Purthermore, in 
examples of the use of this procedure (MSSs 1, 2 
and 3 and SPOT) , the calibration beam passes 
through only a small portion of the aperture of the 
systm, thus not simulating the actual operation of 
the system. When imaging the ground, the system 
entrance aperture is irradiated over its entire 
area by flux incident over a rouchly three 
steradian solid ancle. In the imaging mode there 
is much more stray light present in the system and 
incident on the focal plane. If this additional 
flux level is unknown, it may Introduce a 
substantial uncertainty Into' the absolute 
calibration of the system. 

(c) Reference can be made in-flight to a 
ground area of known spectral radiance. If at the 
time the sensor system is imaging the known area, 
measurements are made of the atmospheric conditions, 
these data can be used with an atmospheric 
radiative transfer program to predict the spectral 
radiance at the entrance pupil of the sensor. The 
main uncertainty in this approach is that of 
determining the atmospheric aerosol content well 
enough. The approach is also limited to scenes 
having large uniform areas of high radiance. Por 
example, although many water bodies are of ' 
sufficient size and uniformity, they are not 
appropriate for calibration purposes because their 
radiance is too low to provide a calibration of 
sufficient accuracy or to cover much of the dynamic 
ranoe of the sensor. Portunately, some suitable 
areas do exist, particularly in the arid reoiona of 



Fig. 3. 


Classification of radioastric calibration procedures. 


173 


3 









the worl^^^ for examole at ''^ite Sands, New Mexico, 
in* the United States. 

The rest of this paoer is devoted to 
discussion of the factor^' and in-fliaht calibration 
of an MLA system usino the self-calibrated 
ohotodiode aporoach and the use of a ground 
reference area for calibration ourposes. 

4. CALIBRATION IN THE FACTORY AND IN ORBIT 

The concept oroposed for the factory 
calibration is similar to the oroposed orbital 
procedure, the main difference being that an 
artificial source is used in the factory and the 
sun is used in orbit— simply a matter of 
convenience in the former case and of convenience 
and reliability in the latter case. In the 
factory, redundancy is not at a premium and our 
requirements for a source are simply power, 
soectral flatness, and stability. We do not need a 
standard source although an array of standard NBS 
FEL tungsten halogen lamps could be used, if their 
polarization characteristics can be tolerated (Ref. 
13) . A xenon arc selected for minimum arc wander 
and with a highly stable power supply and a 
feedback loop would suffice. 

The source would be used to irradiate a 
near-Lambertian, near-unity-reflectance# 
white-surfaced panel perhaps ImxO.Smin size in 
front of the system. (An integrating sphere could 
be used, but it would have to be very large# and 
uniformity checks can sometimes themselves 
introduce non-uniformities.) A self-czilibrating 
NBS-style radiometer# with the incident beam 
perhaps defined by two or three apertures# and 
using spectral bandpass filters matching those used 
in the MZ«A, %fould be used to determine the radiance 
of the panel in each band. The MLA would image the 
panel out of focus# but being an extended object# 
its image would have exactly the same ir radiance in 
or out of focus. The arranaement is sketched in 
Fig. 4. 



Fie. 4. The factory procedure for absolute 
soectroradiometric calibration. 


To avoid problems due to the nonuniform 
irradiation of the panel# the MLA should be rotated 
to sequentially irradiate the focal plane with the 
image of the same small area that is sampled by the 
radiometer. The reason for a large panel is to 
simulate the viewing conditions from apace in 
which# depending on the baffle design# significant 
out-of- the- field-of- view stray light could be 
incident on the image plane to modify the 
calibration. For this reason it would be 
worthwhile to conduct at least one calibration 
using a white panel several meters in diameter. To 


check for linearity of response, several different 
Irradiance levels on the panel should be used. 

This Irradiance level can be changed conveniently 
by inserting heat-resistant neutral density filters 
in front of the stabilized xenon arc source. 

The profiles of the spectral filters should be 
measured in a spectrophotometer using the same 
F/No. beam as the MLA. If they are integral with 
the array, they should be measured using a double 
monochromator, again with the MLA F/No. , before 
installation in the focal plane. Care should be 
taken to cover the whole wavelength sensitivity 
range of the detectors, the off-band suppression 
being particularly important for detectors with the 
wide spectral response of silicon. 

The procedure proposed here for in-flight 
calibration is similar to the panel method sketched 
in Fig. 4, but it uses the sun as the source. Ve 
believe that the irradiance over the panel can then 
be considered to be uniform and known spectrally to 
better than 1% absolute. (Several solar 
measurement programs are currently being conducted 
with this accuracy as a goal. However, if the 
uncertainty is thought to be greater than 1%# a 
pyroelectric detector could be used to measure the 
direct solar flux in orbit# over the wavelength 
intervals of interest and at the same time that the 
system is being calibrated.) The calibration would 
be carried out in the few minutes while the 
spacecraft is sun- illuminated but before it images 
the sun-illuminated earth. The absolute radiometer 
containing the pyroelectric detector would now be 
needed only to check for any deterioration in the 
reflectance of the panel owing to exposure to the 
space environment and short exposures to 
unattenuated 0V and other high energy radiation 
from the sun and from space. In this last respect# 
the panel would usually be stowed in a 
well-shielded compartment and exposed only during 
actual calibration checks. Also# when deployed# it 
would not interfere with normal operation of the 
system# as it would be viewed by the stereo mirror 
in one of its extreme positions. In this respect 
it is fail-safe. 


5. WRITE SANDS AS A CALIBRATION RBFEHBNCE 

To keep the decription of the theory of the 
method brief# spectral dependencies are not 
included in the following discussion* In all cases 
the spectral dependence • is implied# the spectral 
value having the same wavelength dependence as the 
spectral response of each band of the sensor being 
calibrated. The theoretical basis for the 
calibration method is straightforward. 

The radiance at the sensor# L^# is determined 
with respect to the ground radiance# L^# the 
reduction of this radiance by the upward path 
through the atmosphere# sec 9 (where ^*ext 

the atmospheric extinction and 0 is the nadir scan 
angle)# and L^#<^# the atoiospheric path radiance 
(where ^ is the azimuth angle of the observation) • 
Thus# 

^9 ■ ssc6> + 

The path radiance term# Lq#^# In the equation 
is unknown and can be determined onlv by the use of 
a radiative transfer calculation. Such 
calculations require a knowledge of ground 

reflectance# the irradiance incident at the top of 
the atmosphere# the aerosol optical deoth and phase 
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function# the ootical deoth due to otone and water 
•’^apor# the solar zenith angle# the nadir ancle of 
the sensor# and its azimuthal ancle with respect to 
the sun's direction. These radiative transfer 
calculations assume an infinite# flat, Lambertian 
reflecting surface with an atmosphere comoosed of 
plane homogeneous layers. Althouoh these 
assumptions are never exactly met, thev are closely 
aooroximated hv the nearly Lambertian flat surface* 
of the gvosum at White Sands and the small field of 
view of the soace sensors. 

The cuantity t* can he determined from the 
solar radiometer data in the following way. The 
solar radiometer records the irradiance at the 
cround as a function of the solar zenith ancle 
during the morning of the overpaea. If the 
atmosoheric conditions are steady durlnc that 
period# the Langley plot of In against sec6 is 
a straight line of slope r * The extrapolation 
of the line to a aece^ value of zero gives the 
value of the Irradiance above the atmosphere in the 
measurement pasaband. The value of the cround 
radiance Lq is determined at the time of the 
overflight by the solar radiometer looking at the 
ground at the sane B and angles as that of the 
space sensor « The ground radiance value can be 
converted into the reflectance reguired as input to 
the radiative transfer program in one of two ways. 
First# the global irradiance can be measured by the 
solar radiometer with either a horizontea cosine 
diffuser or an integrating sphere# with horizontal 
entrance port# placed over the entrance aperture of* 
the solar radiometer. Or second# the solar 
radiometer can be pointed down at a horizontal# 
diffusing# white surface such as BaSOt* or Halon. 

The cosine diffuser approach is preferred because 
its small surface can be conveniently stored to 
minimize contamination. Intercasparisons will be 
made# during the six monthly recalibrations of the 
field equipment# between the cosine diffuser global 
irradiance and the Lambertian white reference 
surface measurements of the ground reflectance. 

The nadir angle from the satellite to the ground 
observation point# the solar zenith angle# and the 
azimuth of the satellite nadir angle to the sun*s 
direction are all readily determinable. Such 
quantities aa atmospheric pressure# humidity# 
surface and air tesiperature# etc. # are routinely 
monitored by the Atmospheric Sciences Laboratory at 
White Sands. These data will he formatted with the 
radiation data and telemetered to Wallops Island . 
Virginia, united States. 

The difficult quantities to determine for 
radiative transfer calculations are the 
characteristics of the atmospheric aerosols. 
Fortunately# we can make reasonable assumptions 
regarding their values based on prior results. For 
example# 79 aerosol profiles have been measured at 
White Sands (Ref. 14) to establish an atmospheric 
models the effects of aerosol complex refractive 
index and size distribution on extinction and 
absorption for the wavelength range 0.55 to 10.6 urn 
for desert atmospheres are reported in reference 
15? and reference 16 reports that the imaginary 
refractive index is 0.007 at 0.6 urn and shows 
little dependence on wavelencth over the range 0.3 
to 1.1 urn. 


We intend to Investigate sene new techniques 
for determining aerosol chacteri sties that are 
under development and that make use of sky 
polarization measurements. We also Intend to 
explore methods of monitoring the aerosol size 
distribution from inversion of the spectral optical 


deoth measurements (Ref. 17) and aureole 
measurements (Ref. 18) . Some of these 
measurements and analyses# particularly the 
polarimetric# fall into the research category. 

They reoresent attempts to improve the accuracy of 
the calibration procedure and# because of their 
untested nature, recuire validation. However# 
without them we anticipate an overall accuracy of 
about =3% absolute. This represents an improvement 
over the ±5% value claimed in reference 19, and is 
due to our anticipated ±1% absolute accuracy of 
calibration of the ground instrumentation. 

The proposed basic design of the 
Instrumentation is similar to that being designed 
and built at the University of Arizona for the 
investigation of solar irradiance variations over a 
23-year period. It employs a precision alt-azimuth 
tracking stand, with stepper motors to drive the 
two axes# so that it can be pointed in almost any 
direction or be held in alignment with the sun. A 
microprocessor based co m puter system will be used 
to control the motors as well as the data 
acouisition and processing system. A conceptual 
view of the instrument is shown in Pig. 5# further 
details of its design will be presented at the 
symposium. 



Fig. 5. Conceptual view of spectroradioa- 
eter for ground and amospheric measurements • 
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